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Abstract 

We study neutrino oscillation within the framework of three generations in matter. We propose a 
simple method to approximate the coefficients A, B and C which do not depend on the CP phase S in 
the oscillation probability P(u e — > v^) = A cos S + B sin S + C. An advantage of our method is that an 
approximate formula of the coefficients A, B and C in arbitrary matter without the usual first order 
perturbative calculations of the small parameter Am^i/Amli or sin #13 can be derived. Furthermore 
we show that all the approximate formulas for low, intermediate and high energy regions given by 
other authors in constant matter can be easily derived from our formula. It means that our formula is 
applicable over a wide energy region. 

1 Introduction 

Recent experiments clarified that the solar neutrino deficit and the atmospheric neutrino anomaly are 
strong evidences for the neutrino oscillations with three generations. The solar neutrino deficit is explained 
by 24 — > Vy, oscillation pQ and the atmospheric neutrino anomaly is explained by — ► v r oscillation [2]. 
In the recent SNO [3j and KamLAND experiments @], the solar neutrino problem has been solved by the 
large mixing angle (LMA) MSW solution 5 . Furthermore the upper bound of #13 is given by the CHOOZ 
experiment 0. Thus, there are two small parameters 

a = Amlx/Amii ~ 0.03, sin 9 13 < 0.16. (1) 

The remaining problems are the determination of sign Am| 1; the measurement of the 1-3 mixing angle #13 
and the CP phase 5 In the limit of vanishing mixing angle #13 or vanishing mass squared difference 
Am^x, the CP violating effects in the oscillation probability disappear. Therefore the magnitude of the 
two small parameters a and sin 6*13 controls the magnitude of the CP violation. The LMA MSW solution 
in the solar neutrino problem has opened the possibility of the observation of CP violation in the lepton 
sector. For this purpose, many long baseline neutrino experiments are planned [Sj. 

The matter effect received from the earth is important in the long baseline neutrino experiments, because 
fake CP violation is induced due to matter effect. The Preliminary Reference Earth Model (PREM) is 
well known as the model of the earth density and is usually used in analysis of long baseline experiments. 
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However, it has recently been pointed out in geophysical analysis of the matter density profile from J-PARC 
to Beijing [5] that the deviation from the PREM is rather large. In this paper, we derive an approximate 
formula of neutrino oscillation probability without assuming any specific earth density models. 

In constant matter, various approximate formulas have been proposed in low energy 111! |I2 | . in 
intermediate energy ^] and in high energy regions ^3 El EH- In the case that the matter 
density is not constant, approximate formulas have been also derived in ^5] EI3 1^ EH E3 EI3 usm S 
perturbative calculations to analyze the terrestrial matter effect. However, the question of how to separate 
the genuine CP violation due to the leptonic CP phase from the fake CP violation induced by matter effect 
has not been investigated sufficiently in arbitrary matter. 

The next step is to analyze the CP violating effects in more detail in the case of non-constant matter 
density. In order to obtain a hint for this problem, we will briefly review the approach applied in the solar 
neutrino problem. It is difficult to derive the exact solutions for solar neutrino problem in three generations 
except for some special matter profile. As an approach to derive the neutrino oscillation probability, a low 
energy approximate formula was proposed in |26| . By averaging Am^, they derived the formula 

[y e ^u e ) = cos 4 fl 13 P (2) (v e -> u e ) + sin 4 6 13 . (2) 

This is a formula to reduce the calculation of the survival probability P*- 3 -* (y e — > v e ) in three generations to 
that of P( 2 >(v e — > v e ) in two generations. Therefore this formula is called the reduction formula 27 . This 
reduction formula is useful for the analysis of solar neutrino experiments, but it is not directly applicable 
to long baseline neutrino experiments planned in the future, because we cannot average Am^ in long 
baseline experiments. Therefore wc need to derive the reduction formula which is valid without averaging 
Am 2 3l . 

In a series of previous papers we have calculated the oscillation probability P(y e — > v^). In the papers 
|28j we have shown that the CP phase 8 dependence of P{y e — > v^) in constant matter is given in the form 

P{ve -> Vp) = Aeon 5 + BsinS + C (3) 

and have derived an exact but simple expression for the coefficients A, B and C . In the next paper |29j we 
have presented a simple and general formula which does not depend on the matter profile. As a result we 
have concluded that the equation @ is valid even in arbitrary matter. However, in the case of non-constant 
matter density, there exist no closed-form expressions for the coefficients A. B and C . 

In this paper we propose a simple method to derive the approximate formula of the coefficients A, B 
and C taking account of the small parameters a and sin #13 . The coefficients A and B are linear in a and 
sin #13. These coefficients represent the genuine three flavor effect. Therefore it has been considered that 
the first order perturbative calculations of a or $13 are needed for the derivation of A and B. However it is 
possible to calculate A, B and C without the usual first order perturbative calculations of small parameter 
a or sin 0i3 in our method. As we shall see later in section 2, the reduction formula in arbitrary matter is 
derived as 

A ~ 2Re[S%S* e ]c 23 s 23 , (4) 
B ~ -2Im[S^ e ] C23S23 , (5) 

c * \sU?& + W*L, (6) 

where and S!? e are the oscillation amplitudes calculated in the following Hamiltonian, respectively 

H l = O 12 diag(0,A 21 ,A 3 i)Of 2 +diag(a(t),0,0), (7) 
H h = O 13 diag(0,0,A3i)Of 3 +diag(a(i),0,0). (8) 

Here Ay is defined by Ay = Am 2 j/2E, Oij is the rotational matrix in the ij plane and a(t) is the matter 
potential. Since both H e and H h are Hamiltonians in two generations, the equations Q-||HJ| are formulas 
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in order to reduce the calculation of the coefficients A, B and C in three generations to the oscillation 
amplitudes in two generations. Furthermore we show that other approximate formulas known in constant 
matter can be easily derived from our formula. This means that our formula is applicable to a wide energy 
region. 



2 New Idea for an Approximate Formula 

In this section we propose a new idea to derive an approximate formula for neutrino oscillation prob- 
ability. At first we review a general framework for the oscillation probability in arbitrary matter. Next 
we introduce how to derive an approximate formula from this framework. We also discuss the difference 
between our method and usual methods. 



2.1 Review of General Formulation 

In this subsection we briefly review that the CP dependence of the oscillation probability P{v e — > v^) 
is given in the form as P{v e — » v^) = A cos 6 + BsinS + C in arbitrary matter. More detailed calculation 
has been given in the papers 2!) . 

The Hamiltonian in matter is given by 

H = C/diag(0, Aai, Asi)^ + diag(o(t), 0, 0), (9) 
where U is the Maki-Nakagawa-Sakata (MNS) matrix PHj. Using the standard parametrization 

u = o 23 vo 13 r^o 12 , (io) 

the Hamiltonian is written as 

H = O 23 rO 13 O 12 diag(0, A 21 , A 31 )Of 2 Of 3 rto£, + diag(<z(f), 0, 0), (11) 

where T = diag(l, l,e ) is the phase matrix, It should be pointed out that the matter potential of the 
Hamiltonian contains only the e-e component. The Hamiltonian can be written in the form decomposing 
the CP phase 5 and the 2-3 mixing angle # 23 as 

H = 23 TH'^0^, (12) 

where H' is defined as 

H' = O 13 O 12 diag(0,A 21 ,A 31 )Of 2 Of 3 + diag(a(i),0,0). (13) 
The amplitudes S and S' are given by substituting the relations ifT^I and (fT*5f> into the following equations 



S = Texpj-i^ ff(t)cftj , S" = Tcxpj-i^ H'{t)dt\. (14) 

Then we obtain the amplitude for vL — > v' a as the a-(3 component 

S a p = (O^TS'T^O^ais. (15) 
In particular when we choose \x and e as a and /3, the amplitude S^ e is given by 

5 Me =S' c 23 + S' Te s 23 e l5 . (16) 



From this relation, the probability is calculated as 



P(v e -> Vy) = AcosS + Bsind + C, (17) 
A = 2Re[S'; e S' Te }c 23 s 23 , (18) 
B = -2Im[S; e S; e ]c 23S 23, (19) 

J fj,e\ L 23 T" \ u re\ °23i 

which is the exact formula in arbitrary matter derived in the previous paper |29| . 



C = \S' | 2 cl 3 + \S' Te \ 2 sl 3 , (20) 



2.2 Order Counting of A, B and C on a and sin 8i 3 

In this subsection we study how the coefficients A, B and C defined in l)18 [l -l)2() [l depend on a and 
sin 0i3. Instead of A, B and C, we study the dependence of S' and S' Te on a and sin #13 by taking the 
limit either #13 — > or a — * 0. 

At first, taking the limit #13 — * 0, the Hamiltonian reduces to 

H £ = lim Zf' (21) 
= O 12 diag(0,A 21 ,A 3 i)Of 2 +diag(a(t),0,0) (22) 



( 



A 2 isl 2 +a(t) A 2 isi 2 ci 2 
A 2 isi 2 ci 2 A 21 c? 2 J . (28) 
A 3 i 



This Hamiltonian expresses the fact that the third generation is separated from the first and the second 
generations. We simply obtain the amplitude 

S're = (24) 

from the Hamiltonian . It means that the order of S' Te is given by 

S; c = 0(Binfl 13 ). (25) 
for the case 6*13 7^ 0. In the same way, taking the limit A 2 i — > 0, the Hamiltonian reduces to 

H h = lim H' (26) 

= O 13 diag(0,0,A3i)Of 3 + diag( a (t),0,0) (27) 

/ A 31 s 2 13 + a(t) A 31 s 13 ci3 

I . (28) 



\ A 3 isi 3 ci 3 A 3 i< 
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This Hamiltonian expresses the fact that the second generation is separated from the first and the third 
generations. We simply obtain the amplitude 

S; e - (29) 
from the Hamiltonian Q28[). It means that the order of S' is given by 

S; e = 0(a) (30) 
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for the case a 7^ 0. Finally, we conclude that the dependence of the coefficients A, B and C on a and 
sin 6*i3 is given by 

A = 2Re[5^5 , ; e ]c 2 3S23 = 0(asin^ 3 ) ) (31) 
B = -2Imls'; e S' re \c 23 s 23 = 0(asm6 13 ), (32) 
C = \S^\ 2 c 2 23 + \S' Te \ 2 s 2 23 = 0(a 2 ) + 0( S in 2 6 13 ). (33) 

Since both A and B vanish in the two flavor limit, either a — > or sin #13 — > 0, this fact represents the 
genuine three flavor effect. The coefficients A and B are doubly suppressed by these small parameters a 
and sin ^13. In Refs. this is pointed out for the case of constant matter density. However these 

results (|31f) - H33f) are correct even in arbitrary matter profile. 



2.3 Main Result 

In this subsection, we propose a simple method to approximately calculate the amplitudes S'^ e and 
S' Te . From the result of the previous subsection, the dependence of S'^ e and S' Te on a and sin #13 is given 

by 

5^ = 0(a), S; e = O(sin0 13 ). (34) 
We expand both S'^ e and S' Te in terms of two small parameters a and sin 9i 3 as 



S'^ = (o(a) + 0(a 2 ) + 0(a 3 ) + ■ ■ + (0(a sin 6» 13 ) + 0(a 2 sin 6> 13 ) + • ■ ■ ) (35) 

= S^ + Oia sin 9 13 ) + 0(a 2 sin 6> 13 ) + • • • , (36) 

S*; e = (O(sin0i 3 ) + 0(sin 2 6>i 3 ) + •••) + (o(q sin 6>i 3 ) + 0(a 2 sin 9ja) + • • •) (37) 

= 5^ + 0(asin6»i 3 )+0(a 2 sin6l 1 3) + ---, (38) 
where S^ e and S 1 ^, are defined by 

S e ue = lim S' (39) 

S* e = UmS;, (40) 



From (|36fl and l|38[> we can approximate the amplitudes as 

S'lie — ^/j,ei (41) 

S' Te ~ S T fe e . (42) 

The accuracy of this approximation is determined by the magnitude of the higher order terms on sin #13 
and a. At present, the upper bound of sin #13 is given by the CHOOZ experiment. In future experiments, 
when the value of 6*13 will become smaller, the accuracy of the approximate formula can be better. It 
is noted that the simple method introduced in this subsection does not depend on whether the matter 
density is constant or not. We obtain the oscillation probability from the reduced amplitudes as 



P{v e ^v p ) = A cos (5 + £sin£ + C, (43) 

A ~ 2K e [S l ; e S h Te ]c 23 s 23 , (44) 

B ~ -21m[S%S!? e }c 23 S 23 , ( 45 ) 

C * Kf4 3 + \S^s 2 23 . (46) 



5 



This formula is one of the main results obtained in this paper. The advantage of this formula is as follows. 
First, this formula is derived by using only two small parameters a and sin #13 without assuming a specific 
matter density model. Therefore, this formula is applicable to the case of the PREM, akl35f and so on. 
Second, the reduction formula applied to the solar neutrino problem is valid only in the case that the 
averaging for Am^ is possible. However, our formula is effective even in the case that the oscillation 
probability can not be averaged. Namely, it is applicable to long baseline experiments. Third, we derive 
this formula without the first order perturbative calculations of small parameter a or sin #13. This is the 
reason why our derivation is easier than usual perturbative methods proposed by other authors. More 
detailed discussion is given in the next subsection. 

2.4 Comparison with Usual Perturbative Calculations 

In this subsection, we compare our method with usual perturbative methods and describe the advantage 
of our method clearly. From the result of the previous subsection, the dependence of the coefficients A 
and Bona and sin# 13 is given by 

A = 0(asin#i 3 ), B = 0(asin#i 3 ). (47) 

As both a and sin #13 are small parameters, there are two kinds of perturbative methods. One method is 
to consider a as a small parameter and treat #13 exactly. Another method is to consider sin #13 as a small 
parameter and treat a exactly. The former case means that we consider H' as a perturbation from H l 

H' = H l + O(sin0i 3 ). (48) 

We need to perform the first order perturbative calculation to obtain A and B in this perturbative method. 
Similarly, the later case means that we consider H' as a perturbation from H h 

H' = H h + 0(a). (49) 

In order to calculate A and B, we need to perform the first order perturbative calculation. In both cases, 
we need to perform the first order perturbative calculation, because the CP violating effects disappear in 
the limit of vanishing a or #13. 

Let us interpret the above usual perturbative method by using the general formulation ((TBI) and (|T9"|) as 
follows. The expressions A and B are represented by two kinds of amplitudes 5^ e and S' Te . The dependence 
(|25ll and (f of the two amplitudes on a and sin #13 is rewritten as 

S; e = 0(a sin 9 13 ), S' Te = 0(a° ahx9 13 ). (50) 

If we use the expansion in terms of sin #13 to calculate both amplitudes S' and S' Te , the amplitude S' 
can be calculated in the zeroth order perturbation. However, S' Te need to be calculated in the first order 
perturbation. In the same way, if we use the expansion in terms of a, the amplitude S' Te can be calculated 
in the zeroth order perturbation, but £' need to be calculated in the first order perturbation. 

An advantage of our method is that we are able to calculate both the amplitudes S^ e and S' Te in the 
zeroth order perturbation, namely without the first order perturbation of the small parameter a or sin #13. 
If we expand S'^ e in terms of sin #13 instead of a, we do not need to perform the first order perturbation. 
Similarly, if we expand S' Te in terms of a instead of sin #13, we do not need to perform the first order 
perturbation to calculate S' Te . One of the essential points of our method is that the Hamiltonian to 
calculate is different from that to calculate S^: e . Another point is that we only have to calculate S^ e 
and S^,, by using the Hamiltonian H l and H h in the framework of two generations, respectively. These 
ideas make the calculations of the probability easy. 



3 Approximate Formula in Vacuum or in Constant Matter 



In this section, we calculate the concrete expressions for A, B and C both in vacuum and in constant 
matter by using the new method. Moreover, we compare the value of these coefficients with exact value 
by numerical calculation. 

3.1 In Vacuum 

At first, we calculate S^ e in vacuum, namely in the case of a(t) = 0. S^ e is 

SU = (expHff £ L)) Me (51) 
= (0 12 di a g(l,e- <Aali ,e- iA » li )Of 2 ) Me (52) 

= —i sin 26*12 sin — ^ — exp ^— i— ^~L^J (53) 

from the Hamiltonian l|23|) . Similarly S% e is 

S* e = -iam2B la Bw^^exp^-i^Lj (54) 

from the Hamiltonian J2HJ- We obtain the expressions for A, B and C 

A ■ nn ■ ■ A 2 iL , A 31 L A 32 L 

A ~ sin 20i2 sin 2023 sin 201.3 sin — ^ — sin — - — cos — - — , (55) 

A 2 iL A 31 L A 32 L , . 

B ~ sin 20i2 sin 2023 sin 20i3 sin — - — sin — - — sin — - — , (56) 



C ~ c\ 3 sin 2 20 12 sin 2 + s| 3 sin 2 20 i3 sin 2 , (57) 

by substituting and (JS1J into gSJ> an d EHl- 



3.2 In Constant Matter 

Next we calculate the amplitudes in constant matter, namely in the case of a(t) = a. At first, we 
diagonalize the Hamiltonian ill' Ml) in constant matter by the orthogonal matrix 0[ 2 as 

H e = Oi2diag(0,A 2 i,A 3 i)Of 2 +diag(a,0,0) (58) 
= 0^diag(A £ i,A f 2,A 31 )(0^) T (59) 
to calculate Sw Here \\ (i = 1,2) is the eigenvalue given by 



Af = | ( A 2 i+o± V /(A 2J «,,s20 I2 -■■ + A], sin 2 20, 2 j . (OOi 

and Af and A2 correspond to the sign — and the opposite sign +, respectively. The effective mixing angle 
sin 20i 2 is calculated as 

sin 2 20^2 = Al^26 12 

(A 2 i cos 20i2 - a) 2 + A 2 ! sin 2 20 i2 
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From iHiOll and l|61(l we obtain the relation 



A e 21 _ sin26»i2 / / a 



2 



cos 2012 - ) + sin 2 29 12 . (62) 



A 21 sin 26^f 2 V V A 2i 
The amplitude S^ e is calculated by using the Af and sin 20f 2 as 



Si = (exp(-i^i)) Me (63) 
= (o 12 diag(e- lA i i , e -^ i ,e- iA3li )Of 2 N ) (64) 



- i sin2^ 2 sin^|^exp(- i ^i^ 2 -Lj (65) 



-z sin 20£ 2 sm — Tp - ex P I ~ l 2 ' ' ^ ' 



Next let us calculate S^ e from the Hamiltonian (|28|l diagonalized by the orthogonal matrix Of. 



13 

= O 13 diag(0,0,A 3 i)Of 3 + diag(a,0,0) (67) 
= O^ 3 diag(Aj,0,A5)(O^) T . (68) 

The eigenvalue = 1, 3) of this Hamiltonian is given by 

Af = i ^A 3 i + a ± ^ (A 3i cos 20 i3 - a) 2 + sin 2 20^ , (69) 

where Af and A3 correspond to the sign — and the opposite sign +. Moreover we obtain the effective 
mixing angle sin20f 3 is calculated as 



sin 2 2** = A ^ lSin 2dl \ , . (70) 

(A 3 iC0s26»i 3 -a) 2 + A 2 1 sin 2 26» 13 



From (|69|) and 1)70(1 we obtain the relation 



AS, sin26»i 3 // a A 2 , 

- -''cos 20i3- -r— + sin 2 20i3. (71) 



A31 sin 20^ V V A 3i 
The amplitude S^, is calculated as 



^ e = -i sin 20^ sin ^ exp f -i ^tf^ . (72) 



Substituting and (O into (JSJ), gSJ and @BJ), we obtain 



A ~ sin 26\ 2 sin 20 2 3 sin 20^ sin sin cos — ) 



(73) 
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B 



C 



A 1 L A h 

sin 20i2 sin 29 23 sin 29\ 3 sin — ^— sin — 



L . A 32 L 

sin , 

2 2 ' 



cf 3 sin 2 20f 2 sin 2 ^ 



sin 2 26# 
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sin' 



A h 31 L 



(74) 
(75) 



The low and high energy MSW effects are contained in sin 29\ 2l ^21 an d m sm 2^13, A^ of the approximate 
formula, respectively. This is the reason why this approximate formula is applicable to a wide energy region. 
The term including B, which is proportional to sin 5, is related to T violation |S1 E2 El ■ However 
it is difficult to observe only this term in future long baseline experiments. Therefore there are many 
attempts to extract the information on the CP phase from the terms including both the coefficients A and 

b EU EU ET|. 

Next let us compare our approximate formula with the exact one. We use the parameters Am 21 = 
7.0 x 1(T 5 eV 2 , Am 2 31 = 2.0 x 10~ 3 eV 2 , sin 2 29 12 = 0.8, sin 2 29 23 = 1, sin#i 3 = 0.16, the oscillation 
length is L — 730 km a = y/2GpN e , where Gf is the Fermi constant and N e is the electron density 
in matter calculated from the matter density p = 3g/cm 3 and the electron fraction Y e — 0.5. We plot 
the coefficients A, B and C as a function of the energy within the region O.Of GeV < E < lGeV. These 
coefficients calculated from the exact formula are compared with those from the approximate formula in 
Fig. 1. From this figure we find that the approximate formula almost coincide with the exact formula. The 
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Figure 1: Comparison of our reduction formula with the exact one in the coefficients A, B and C, plotted 
from top to bottom and the exact, approximate formula and the difference from left to right. 



error is estimated to be less than 20 % from Fig. 1. The difference between the exact and the approximate 
formulas is caused by ignoring the higher order terms in the perturbative expansion on sin $13 and a. 
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4 Derivation of Other Approximate Formulas 



In this section, we derive the approximate formulas given by other authors in constant matter from 
our formula. There are formulas for the low energy |X(JI 111! IT2] , the intermediate energy \l'6l 1141 ITS] and 
the high energy regions ^3 El EH • 

4.1 Low Energy Formula 

At first we derive a low energy formula with large mixing angle #12, which is similar to those in 
[1011111 IT2*] . Under the low energy condition 

a«A 3 i, (76) 

the following relation 



Ah 
A31 



sin 26*13 
sin 20f 3 



(77) 



is derived by expanding A 31 and sin20^ 3 in terms of a/A 31 . Namely A 31 and sin 20^ 3 in matter can be 
approximated by the quantities in vacuum. Furthermore, if we take the limit O12 — > 7r/4 



A|i _ sin26>i 2 _ yg| + 
A 2 i ~~ sin20f 2 ~ A 2 i 



(78) 



is obtained. By using the relations (f77|) and (JTSJ, the coefficients A,B and C for (|75|) . l (74*j l and ((ToTl are 

reduced to the following expressions 

a A 2i -oa -on -on • . + a 2 L A 32 L 
' sm 20i2 sm 2023 sm 20i3 sm — 7: — sin- — £i - cos — - — , (79) 



2 io ±j 2 2 2 



B ~ - =sin26>i 2 sm 2023 sm 20i3 sm sin- sin— — , (80) 

a/A^ + a 2 2 2 2 



C 



A 2 
^21 



A 2 
^21 



23 



sin 2 20i2 sin' 



2 L 



s 23 sin 20-| 3 sin 



A31L 



13 ( 



(81) 



where the condition derived from the low energy condition i|76|l 



sm — - — ~ sin (82) 

2 2 y J 



is also used. The applicable region for energy is given by 



Am 2 i \ /3g/cm 3 



E I:,(>v f l ~ — 1 

from the condition 176(1 . In addition to this condition, the applicable region of (|79|l . l(80[l and (jSljl is 

restricted by 

1<<8000tm (_|_)(l£l£X!V m 
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which comes from the approximation in the oscillation parts of l|79[l. (|80|l and l|81|l. A similar result can be 
obtained for the perturbation of sin0 13 |19) . They have proposed a low energy formula in arbitrary matter, 
by using the first order perturvative calculations. Our method has the advantage that the calculation is 
much simpler. This approximate formula coincides with that in vacuum in the low energy limit, or in other 
words, this result recovers the vacuum mimicking phenomenon which has been discussed in |12l I38| . 

4.2 Intermediate Energy Formula 

At first we derive the intermediate energy formula jlHI 1141 115) from our formula. Under the low energy 
condition 



a< A 3 i, 

we expand A 31 and sin20^ 3 up to first order of a/A 3 i 

A 31 ~ A 3i - 2a cos 26»i3, 

2a 



sin26> 



13 



sin 26»i3 ( 1 + -7 — cos 26>i 3 
A31 



(85) 

(86) 
(87) 



Substituting QJgft and (JH3 into (O, we obtain the expression 



15: 



h |2 
re I 



s 2 23 sin 2 26»i3 ( 1 



2a 
A3T 



cos 26*13 I sin 



(A31 - 2acos26»i 3 )L 



S23 sin 2 26* 



13 



-j^— cos 26*13 ] sin 2 ^77— — aL cos 26*13 sin(A 3 iL) 
A31 / 2 



(89) 



where we also use the approximation 



aL cos 2^13 -C 1 



(90) 



from the first line to the second line. Furthermore, under the assumption that A21-L/2 is small, we can 
approximate 



ALL ALL A 2 iL 

sin ^ Hr" ^ 

2 2 2 

and from (|73)l . (|74(l . (|75(l . the approximate formula for A, B and C is derived as 



(91) 



.4 



B 



C 



1 A 2 iL 

- sin 26»i2 sin 26» 2 3 sin 26»i 3 - sin( A31L) , 



• nn ■ nn ■ no ^2\L . 2 A 31 L 

sm 20i2 sm 2023 sm 2t/i3 — ^ — sin — - — , 



s 2 3 sin 2 26»i3 



1 + — — cos 26*13 ) sin 2 — — aL cos 26*13 sin(A3iL) 

A31 / 2 



(92) 
(93) 
(94) 



One of the conditions for the applicable region of this approximate formula, namely for the upper limit 



E < 15 GeV 



Am§i \ {3 g/cm 
IO- 3 eV 2 



(95) 
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is the same as that in the low energy region (|76|) . In addition to this, the conditions due to Q90(l and 1|) 



L « 1700knW 3g/Cm3 ) , (96) 



E » 0.185 GeV ' 21 , (97) 
VlO" 4 eV 2 J V730k m y v ' 

should be satisfied. This approximate formula has been derived by using the perturbations of a and 
a/Amjj |13|. From (|95p. i|96|) and (|97|l . the applicable region is rather restricted because of the expansion 
of the oscillation part. On the other hand, it has the advantage that the contribution of the genuine CP 
violation can be easily distinguished from that of the fake CP violation. 

4.3 High Energy Formula 

Next, we derive the high energy formulas [1611171 IT%| from our formula. Under the high energy condition 

a>A 2 i, (98) 

we obtain 



A 21 sin 2#i2 a 
A 2 i ~ sin26>f 2 ~~ A 2 i 



(99) 



by expanding A 2 i and sin W\ 2 up to the first order of A2i/a. In addition, using the approximation #13 — > 0, 
we obtain another relation 

^L^,! ». (100) 
A 3 i sin 2^3 A 3 i 

The concrete expressions of A, B and C are derived as 

A21A31 . aL . (A 31 -a)L A 32 L 

A ~ — — sm 2t^i2 sm 26'23 sm 26'i3 sm — sin cos , (101) 

a(A3i — a) 2 2 2 

A21A31 . aL . (A 31 -a)L . A 32 L 

B ~ — — sm 2^12 sm 26*23 sm 26*13 sm — sm ^ sin — - — , (102) 

a(A3i — a) 2 2 2 

C ~ gfkcg, sin 2 29 12 sin 2 ± + ^ s 2 23 sin 2 20 13 sin 2 (Asi " " )L , (103) 

2 (A31 — a)- 2 2 

by substituting 199JI and (jl(K)fl into l)73|l . I|74|l and l|75|) . where we also use the approximation 

sin ^"- sm T' (104) 

The applicable region of this approximate formula is calculated from l)98|) as 

<«*> 
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In addition to this, the applicable region is also restricted by 

L<< 8000 km (^)(i^), (106) 

which is derived from l|104|l . Although these high energy approximate formulas (|101|l . H102|) and ()103() 
have been derived at first in |17l IT8*| , their derivation is complicated because of the calculation up to the 
first order perturbation of a. Here, we have presented the simple derivation of these formulas by using the 
new idea of taking only the zeroth order perturbation. 

5 Summary 

In this paper, we study the oscillation probability in matter within the framework of three generations. 
The results are as follows. 

1. We have proposed a simple method to approximate the oscillation probability in arbitrary mat- 
ter. Our method provide an approximate formula in arbitrary matter without the usual first order 
perturbative calculations of the small parameter Am^/Am^ or sin #13. 

2. The concrete expressions for our approximate formula in constant matter has been derived to investi- 
gate the accuracy of the reduction formula l|43l) - (|46[l . We have shown that our formula is numerically 
in good agreement with the exact solution with reasonable accuracy. 

3. We have shown that both the low energy the intermediate energy pU^J^] and the 
high energy |16l 1171 118j approximate formulas in constant matter presented by other authors can be 
easily derived from our formula. This means that our formula is applicable to a wide energy region. 
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